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C\J I There is strong evidence that the mass in the Universe is dominated by dark 

^ ' matter, which exerts gravitational attraction but whose exact nature is un- 

I ■ known. In particular, all galaxies are believed to be embedded in massive 

^ ■ haloes of dark matter. ^'^ This view has recently been challenged by surpris- 

O ■ ingly low random stellar velocities in the outskirts of ordinary elliptical galax- 

r-| ! ies, which were interpreted as indicating a lack of dark matter. ^'^ Here we show 

that the low velocities are in fact compatible with galaxy formation in dark- 
O I matter haloes. Using numerical simulations of disc-galaxy mergers, ^'^ we find 

I that the stellar orbits in the outer regions of the resulting ellipticals are very 

■ elongated. These stars were torn by tidal forces from their original galaxies 

■ during the first close passage and put on outgoing trajectories. The elongated 
^ ■ orbits, combined with the steeply falling density profile of the observed trac- 

. ers, explain the observed low velocities even in the presence of large amounts 

of dark matter. Projection eff"ects when viewing a triaxial elliptical can lead to 
even lower observed velocities along certain lines of sight. 



Introduction 

The common spiral galaxies are known to reside in extended dark-matter (DM) haloes. 
The rotational speeds of their gas discs do not decline outside the visible body,^ unlike 
the expectation from Keplerian circular velocities at a radius r about a mass M, V"^ = 
GM/r. Thus, the DM mass within r is growing roughly as M(r) oc r and it dominates 
the gravitational potential beyond a certain radius. An extrapolation based on the typical 
halo density profile^ found in simulations of the standard ACDM cosmology predicts an 
outer "virial" radius i?vir that is 50-100 times larger than the characteristic stellar radius, 
enclosing 10-20 times more DM than luminous matter, now also indicated observationally..* 
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The conventional wisdom is that the potential wells created by the DM are crucial for 
seeding the formation of galaxies. 

The standard hypothesis is that ellipticals originate from mergers of discs^^ and should 
therefore be embedded in similar DM haloes. There is evidence for DM in giant ellipticals, 
from X-rays^^ and gravitational lensing^^ and in nearby dwarf galaxies from stellar kine- 
matics. -""^ However, ordinary ellipticals lack obvious velocity tracers at the large projected 
radii Vp where the DM is expected to be important. This is typically beyond Rcs,^^ the "ef- 
fective" radius encompassing half the projected light, while measurements of the projected 
velocity dispersion dp of the stellar light are limited to Vp < 2i?efr- 

The strong [OIII] emission line at 5007A from Planetary Nebulae (PN) — hot shells of 
gas expelled from dying stars of (1 — S)Mq — provide a unique tool for extracting crp(rp) 
beyond it!eff- Romanowsky et al.^ measured Up from ~100 bright PNs in each of three 
normal ellipticals, NGC 821, 3379 and 4494, adding to 531 PNs in NGC 4697 by Mendez et 
al..^ They find that cxp typically drops by a factor ~ 1.6 between Vp = R^s and 3-Rcfr- The 
kinematic modelling by the observers^ yielded low mass-to-light ratios, e.g. M/L ~ 7 at 
5i?eff for NGC 3379, consistent with a "naked" stellar population. They interpreted this as 
"little if any DM in these galaxies' haloes" , in conflict with the conventional picture. While 
noticing that increasing velocity anisotropies could in principle produce declining ap, they 
ruled out such "pathological" orbit structure. Similar conclusions were obtained later from 
other ellipticals.^^ The apparent challenge to theory has already triggered radical explana- 
tions.^'^ However, the earlier analysis*^ might have missed alternative solutions because it 
was limited to specific density-profile shapes in stationary spherical systems, the halo PNs 
were identified with the central stellar population, and their maximum-likelihood method 
may suffer from an incomplete orbit library or questionable convergence properties. 

Anisotropic velocity dispersion 

For given density profiles, a lower ap can result from more radial velocities. The dynamics 
implies a lower 3D velocity dispersion a because the pressure needed for balancing gravity 
is provided by a radial (Tj. that corresponds to a lower a. The projection introduces a 
further decrease in dp. This can be illustrated by toy models made purely of circular or 
radial orbits, with the same constant speed and random orientations. If the stellar-density 
profile is steep enough, ap is dominated by the tangential contribution near the equatorial 
plane perpendicular to the line of sight, which is high for circular orbits and low for radial 
orbits. 

The 3D profiles of any component of a spherical gravitating system in equilibrium obey 
the Jeans equation^^ (derived from the coUisionless Boltzmann equation) , 

V'{r) = [a{r)+^{r)-2P{r)]a',{r) , (1) 

a manifestation of local hydrostatic balance between the inward pull of gravity (left) and the 
outward push of pressure (right). Here V'^{r) — GM{r)/r is the squared circular velocity. 
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The steUar density profile i/(r) enters via a = — dlniz/dlnr. Its velocity dispersion consists 
of radial and tangential components, a'^ — + 2(7^; we define 7 = — dlncr^/dlnr. The 
velocity anisotropy is (3 = 1 — ag/a^, with /3 — — oo,0, 1 for circular, isotropic and radial 
orbits respectively. The projection can be performed analytically when /3, a and 7 are 
constant with r (power-law profiles, V'^ — V^^(r/i?eff )"'*'): 

{A in Methods). Note that is a decreasing function of /3 and of a (for a + 7 > 3 
and P > 0) respectively. Local fits to the de-projection of the standard de Vaucouleurs^° 
surface-brightness profile of ellipticals give a ~ 3.13-3.37 at 2-3it!efT-^'^ Our fits to o"p(rp) in 
the observed ellipticals beyond Rcs yield 7 ~ 0.8 ± 0.2. These give A{a,'-f) ~ 0.2 and Cp 
drops by a factor ~ 1.5 between /? = and 1. We learn that one could match the low cxp 
at large radii either by a low Vq or by a high there, and that a high a helps. 

In a more realistic model, we assume a Sersic stellar density profile^^'^^'^^ (= de Vau- 
couleurs' for Sersic index m = 4), a standard DM density profile^ with a typical concen- 
tration (~ 10),^^ and a virial stellar mass fraction ~ 0.06. With P — we recover the 
apparent 5D discrepancy, but m ~ 4 and /9(r > -Refr) — 0.5 (independent of its behaviour 
well inside i?efr) yield a ~ la agreement with the observed dp. A good fit is obtained with 
either m ~ 4 and f3 ~ 0.75, or with m ^ 2.3 (for which a ~ 3.5 near 2.5i?cfr) and j3 ^ 0.5. 
The required f3 is higher than the (3 < 0.2 predicted for DM particles in typical haloes at 
a few i?eff -Rvir)/^ SO the PNs must not trace the DM kinematics. 



Merger simulations 

Assuming that ellipticals form by mergers, and that major mergers of discs can reveal 
generic features of mergers in the ACDM cosmology, we appeal to a suite of simulations 
of such events. ^'^ Two spirals arc put on a parabolic orbit, each consisting of stellar and 
gaseous discs and a bulge, all embedded in a ACDM halo, constructed to match a range of 
typical disc galaxies. The gravitational and hydrodynamical evolution is followed using an 
SPH code,^^ including gas coohng, star formation and supernova feedback (Methods). 

Fig. 1 shows the stacked 3D profiles of the merger remnants. The DM density profile 
is shghtly flatter than an isothermal sphere, p oc r~^, similar to simulated ACDM haloes 
after they have responded to gas dissipation.^^ The robust stellar density falls off more 
steeply, p oc r~^-^, as in elliptical galaxies obeying the de Vaucouleurs profile near 2-3i?eff, 
and with Res ~ 0.015i?vir- For the "young" stars it is somewhat steeper, p oc r~^-^. The 
total-to-stellar mass ratio rises from ~ 2 at 3i?cfr to ^ 14 at i?vir, corresponding at 5i?cfr 
to M/L ~ 15 (compared to the earlier"^ M/L ~ 7) (both for stellar M^^/L = 6). 5D The 
3D a profiles of the DM and stars have similar slopes [as in eq. (2)], falhng off roughly as 
(7 oc r~°-^. 



Our main new finding is the high (3 of the stellar halo velocities. While the DM velocities 
are almost isotropic {P ~ 0.1), the typical stellar P grows from small values at r < Res 
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Figure 1. Three-dimensional profiles of the simulated merger remnants. Ten galaxies at 
two different times after the merger (typically 0.8 and 1.3 Gyr) are stacked together. Shown 
are the profiles for the dark matter (blue) and the stars (red), divided into the old ones 
from the progenitors (dotted) and those newly formed during the merger (dashed). The 
scaling is such that the curves for the stars (all, solid red) are matched at Res- The shaded 
areas mark Icr scatter. The panels refer to density p, mass M and circular velocity V, 
velocity dispersion a and anisotropy /3, with "efF' referring to the quantities at -Refr- 
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(sometimes negative, reflecting a small disc that is irrelevant to the profiles beyond i?cfr) 
to /3 ~ 0.5 at r > Res- In one case P > 0.75, but in another it remains < 0.2. Given V{r) 
in the Jeans equation, the higher a is compensated by a higher (3 and a lower a. 

The simulations demonstrate that the stellar halo originates from tidal processes during 
the first pericentre passage. Some of the halo stars are associated with the two cores and 
the tidal bridge between them; they pass near the center at the coalesence before flying 
outward on radial orbits. Other halo stars first flow out in extended tidal tails and later fall 
back on radial orbits [Supplementary Information, hereafter SI]. Indeed, we find that f3 is 
correlated with the strength of the tidal interaction; it is higher for more head-on collisions 
and when the spins are aligned with the orbit. 

Projected profiles 

The systems are "observed" from three orthogonal directions and stacked together, pro- 
viding a robust average profile and the scatter about it. The data are scaled similarly 
(Methods). Fig. 2 shows the simulated surface-density profile and those of NGC 821,^''' 
3379^^ and 4697.^^ They are all fit by S oc r^^-^ in the range l-5i?efr (as in dc Vaucoulcurs 
profile at ~ 2/?eff)- The simulated projected axial ratios near i?efr range from 1:1 to 1:2, 
and the ellipticity is supported by an anisotropic, triaxial velocity dispersion rather than by 
rotation, similar to ellipticals [SI]. In addition. The distribution of global properties of the 
simulated remnants, such as luminosity, radius and velocity dispersion, is consistent with 
the "fundamental plane" of ellipticals [SI]. Thus, the merger remnants seem to resemble 
typical ellipticals near ~ i?eff in every relevant respect. 

The velocity dispersions in Fig. 2 demonstrate our punchline. While the DM fjp indeed 
lies above the outer observed points, the stellar cTp, ~ 30% lower, provides a good fit 
everywhere in the range 0.5-4i?eff- The slope of the simulated (Jp{rp) in the range 0.5-6i?eff 
is 7 = 0.53 ± 0.16 for "all" and 7 = 0.61 ± 0.22 for the "young" stars, both consistent with 
the 7 = 0.59 ± 0.13 observed for PNs. No fiddling with model parameters is involved — 
simply stacking a sample of merger remnants as simulated. 

The (7p of "young" stars is lower by ~ 9% at 3i?eff (due to their larger a) . Stellar theory 
indicates that these objects, < 3 Gyr old, may represent the observed PNs. Emerging 
from 1.4-2.5M0 stars, the PNs are expected to be much more luminous than those of 
the older, less massive stars, which fall below the detection limits^'^ (Methods). However, 
the radial orbits are a generic result independent of the degree of gas dissipation during 
the merger: our mergers with initial gas-to-baryon ratio ranging from to 70% show 
negligible differences in (3 beyond it!eff- Dissipation results in a more centrally concentrated 
stellar distribution, associated with ^ 10% reduction in outer Up [SI]. We also find that 
the radial orbits and low dp emerge from major and minor mergers alike, independent of 
the progenitor mass ratio [SI], and that the presence of a ~ 22% bulge does not make a 
significant difference. The tidal origin of the stellar halo explains this robustness to many 
merger characteristics. Furthermore, one simulation continued till 3.5 Gyr after the merger 
with no sign of evolution in /3(r) and crp(r) beyond ~ R^s [SI]. 
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Figure 2. Projected profiles: simulated galaxies versus observations. Top: surface density. 
Bottom: velocity dispersion. The merger remnants are viewed from three orthogonal di- 
rections and the 60 profiles are stacked such that the stellar curves ("all") match at Res- 
Colors and line types are as in Fig. 1. The < 3 Gyr "young" stars may mimic the observed 
PNs. The la scatter is marked by shaded or hashed areas or a thick bar. The galaxies 
are marked green (821), violet (3379), brown (4494) and blue (4697) with la errors; PNs 
(circles) and stars (crosses). The surface densities shown for 3 galaxies almost coincide with 
the simulated profile. Green lines refer to the earlier models^ with (upper) and without 
(lower) DM. 
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The ±20% scatter in Up is partly due to the angular momentum of the merger orbit and 
the relative spin inclinations, but also due to the line-of-sight relative to the principal axes 
of the triaxial velocity-dispersion tensor (or rotation axis, A. Burkert et al., in preparation). 
When viewed "face-on" , some of the remnants show cTp lower than observed, while other 
extreme "edge-on" cases show (7p almost as high as that of the DM [SI] . 

The simulated line-of-sight velocity distribution (LOSVD) is consistent with the data 
also beyond the second moment. Both the simulated and observed deviations from a 
Gaussian distribution are typically small, with the fourth moment /14 = 0.03 ±0.05 for the 
central stars [SI]. At r ^ Res, radial, prograde mergers produce negligible values as in 
the PNs of NGC 3379^ (or small positive values as in NGC 5128^^), while more circular 
and retrograde mergers, or gas-rich mergers, can produce negative /i4 as in NGC 4697^ [SI]. 

Conclusion 

We conclude that the PN data are consistent with the simple picture where normal ellip- 
ticals also live in massive DM haloes. The low CTp is primarily due to the radial orbits of 
the halo stars, being tidally ejected from the inner regions during mergers independently 
of dissipation and mass ratio. This generic origin of the radial orbits indicates that the 
results based on our sample of simulations are representative of a broader range of merger 
types expected in the ACDM cosmology, and argues that the low observed PN velocities 
are a natural outcome of this standard model. This should be confirmed by cosmological 
simulations (e.g. low- resolution results,^^ and work in preparation by J. Navarro et al.). 

The range of merger properties leads to a variety of /3 and Up profiles, and the triaxiality 
adds directional variations, allowing extreme (7p values smaller and larger than the PN data. 
The possible association of the PNs with the younger stars, whose density profile is slightly 
steeper, may help reducing their o"p a bit further. Other tracers, involving old stars, are 
expected to show a somewhat higher cTp. This is especially true for globular clusters, ^^'^^ 
given their fiatter density profile^^ and presumably lower anisotropy (due to tidal disruption 
in radial orbits) . A somewhat higher cTp may also be expected in gas-poor elliptical-elliptical 
mergers, common especially in groups (as observed^^), where the collision orbits may be 
circularized by two-body relaxation and dynamical friction. 

Recall^ that a nearly isotropic "naked" stellar system also provides a fit to the crp(rp) 
data inside a few Rcs, appealing to a low Vq rather than a high f3 in eq. (2). Our simulations 
provide another solution, which does include a standard DM halo. The DM model predicts 
that CTp fiattens toward ~ 10i?eflf (as in NGC 5128^^), where the "naked" model predicts 
a continuing decline beyond SRes and very low (7p for other tracers as well. While the 
"naked" model violates much of what we know about galaxy formation and cosmology, the 
DM model, with the radial stellar orbits, seems to be a straightforward outcome of the 
self-consistent picture of structure formation in the Universe. 
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Methods 
Equation 2 

The coefficient in eq. (2) is 

Afr. = 1 r[(a + 7-l)/2] T[a/2] 
^"'^^ (a + 7) r[(a + 7)/2] r[(a - l)/2] " 

Tfiis is a weak function of a and 7. For instance: A{3.5, 1.0) ~ 0.18, A{3.0, 1.0) ~ 0.20, 
A(3.5,0.4) ~ 0.23, A(3.0,0.4) ~ 0.26. 

Merger simulations 

The merger simulations^'^ represent some of the major collisions that likely occurred during the 
hierarchical structure formation according to the ACDM cosmology. The evolution is followed 
using the entropy-conserving, gravitating. Smoothed Particle Hydrodynamics (SPH) code GAD- 
Qgrp 25 Qg^g cooling, star formation and supernova feedback are treated using simplified recipes 
that were calibrated to match observed star-formation rates. 

The progenitor disc galaxies mimic typical big spirals galaxies: one type (G) representing 
today's Sb galaxies and another type containing more gas as in Sbc-Sc galaxies and at high 
redshift. The sample consists of four G mergers, with DM masses 1.2 x (except one 

5 X IO^^Mq), and five Sbc mergers plus one Sc merger with 8 x 10^"^ Mq haloes. The baryonic 
fraction is ~ 5% of the DM halo mass in the G cases, and ~ 13% in the Sbc-Sc cases. The fraction 
of baryons in gas is ~ 20% in G, 52% in Sbc, and 70% in Sc. The particle mass is < IQ^Mq for 
gas and stars and ^ lO^M© for DM. The smoothing lengths are h = 100 and 400 pc with the 
force becoming Newtonian at > 2.3h. 

Two identical galaxies are set on parabolic orbits and merge because of dynamical friction due 
to their massive haloes. Our sample consists of several different orbits and orientations, including 
prograde and retrograde configurations in which the galaxy spins are aligned or antialigned with 
the collision orbital angular momentum. The merger results in two succesive starbursts, one after 
the first close approach, and the other after the second, final coalescence [SI]. The starbursts 
typically occur 1-2 Gyr after the beginning of the simulation, and the remnant is "observed" ~ 1 
Gyr later. The amount of stars formed during the merger is roughly proportional to the initial 
gas fraction, and is not too sensitive to the orbit or orientation. The instantaneous rate varies 
in the range 10-100 Mq/jv. The young stars formed during the merger constitute ~ 30% of the 
total stars; typically 20% in the G cases and 40% in Sbc. The remnant galaxies resemble normal 
elliptical galaxies, as demonstrated above [SI]. 

The smaple used from the suite of 1:1 merger simulations^'^ consists of the following runs: 

The G cases are: G3G3b-ul (the fiducial G case, a gas-poor galaxy obeying the scaling 
relations from SDSS), G3blv5G3blv5b-ul (bulgeless), G3G3r-ul (a retrograde orbit), G2G2r-ul 
(a slightly smaller gas-poor galaxy obeying the SDSS scalings, on a retrograde orbit). 

The Sbc-Sc cases are: Sbc201a-u4 (the fiducial Sbc case, gas rich, small bulge, same as run 
n21ow^), Sbc201-uu43 (with stronger feedback, n2high6), Sbc203-u4 (retrograde orbit), Sbc204-u4 
(very radial orbit), Sbc205-u4 (more tangential orbit), Sc201-u4 (mimicking Sc galaxies with no 
bulge and a higher gas fraction). 
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Scaling the data 

The observed galaxies are presented together using the i?eff determined from each surface- 
brightness profile.28.29,4 We note that i?eff for NGC SSTQ^^'^^ is ~ 50% larger than that quoted.^ 
In Fig. 2, an open circle marks the last point had we used the smaller^ R^s; this uncertainty does 
not make a qualitative difference. The amplitudes of S and dp are scaled by least-squares fits 
of the stellar data at r > 0.2i?eff to the stacked simulated profile as a reference. Replacing this 
reference by a different function of a similar general shape yields similar results. Using only the 
stars at larger radii (up to r > Res), or using the PNs alone, yield Up adjusting factors that differ 
only by a few percent. 

The i?efr of NGC 821, 3379, 4494, 4697 match the mean simulation value after multiplication 
by 0.667, 1.57, 1.00, 1.13, indicating that the simulated and observed galaxies are of comparable 
sizes. The cxp were adjusted by factors 1.00, 1.19, 1.21, 1.11 for best fit. Being comparable to the 
radius scaling factors indicates that the observed and simulated galaxies have a similar velocity 
structure. The mean and la scatter in the simulated remnants are i?eff = 4.05 it 1.04 kpc and 
crp(Eeff) = 154 ± 33kms-^ 

Age of Planetary Nebulae 

The [OIII] luminosity of a PN with mass < 2.5M0 is strongly increasing with the parent stellar 
mass^^ [Figs. 10, 14], hence sharply decreasing with age. A limiting magnitude M5007 then 
corresponds to a maximum stellar age t. For the complete sample of 328 PNs in NGC 4697^ it is 
M5007 ^ -2.6, namely tMoi < 3 Gyr^o [Figs. 18, 19, 26]. With only ~ 100 PNs per galaxy,^ the 
magnitude limit is brighter (by -0.8 magnitudes based on telescope gathering areas), so the stars 
are even younger. Based on theoretical PN luminosity functions^" [Figs. 18, 26], if the population 
is typically older than 1 Gyr then t^^os < 2 Gyr. We therefore adopt t < 3 Gyr as a limit for 
most of the observed PNs in the four galaxies. This indicates an association with the "young" 
simulated stars, and that the mergers of gaseous discs are relevant to those ellipticals showing 
PNs. A caveat is the apparent relative invariance of the PN luminosity function between galaxies, 
seemingly independent of signs for a recent major merger. When there are no signs for a recent 
major merger, the observed PNs, if indeed young, may be the signature of recent minor mergers, 
which are expected to produce similar effects. 
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SUPPLEMENTARY INFORMATION 

The following is an extension of the Letter to Nature, aimed at providing more detailed support 
to the results reported in the Letter. 

Outline: 

L The generic origin of a stellar halo with radial orbits 

2. Robustness to dissipation 

3. Robustness to major versus minor mergers 

4. Robustness to the presence of a bulge 

5. Long-term stability 

6. Parametrizing the velocity-dispersion profiles 

7. The full LOSVD and the /i4 moment 

8. Matching other properties of ellipticals 

9. Scatter induced by triaxiality 

10. On the PN age 

11. Rising versus constant /3 

12. Why was the DM solution missed before? 

13. Robustness and Uniqueness of the results 

Fig. 3 shows snapshots from one of our merger simulations, brought here for general impression 
of some of the issues discussed in the Letter and in this SI. 

1 The generic origin of a stelleir halo with radial orbits 

Fig. 4 demonstrates the important role of tidal processes in the formation of a stellar halo and its 
radial orbits. The final halo stars are composed of two populations: (a) stars that passed through 
the center near the second-pericentre coalescence, some of which were part of the tidal bridge 
formed after the first-pericentre passage, and (b) stars that first flew outward in extended tidal 
tails after the first-pericentre passage and then fell back after the merger of the cores. 

The tidal origin of the radial orbits explains the robustness to the degree of dissipation, the 
progenitor mass ratio, the bulge fraction and so on. It predicts that /3 should predominantly be 
a function of the strength of the tidal interactions associated with the collision. Indeed, we find 
in our simulations a strong tendency for high (3 in collisions of low impact parameter, as well as 
in prograde collisions, where the tidal effects arc stronger. High angular momentum collisions, 
and retrograde collisions, typically lead to low beta values. Such collisions may be responsible for 
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Figure 3. Visual impression. Stellar density maps in one of the merger simulations (fiducial 
Sbc). "All" the stars (blue) and the "young" stars formed during the merger (yellow 
contours). The time sequence (time marked in Gyr) is viewed "face-on" in the plane 
defined by one of the initial discs and the two centers of mass. Young stars form especially 
after each of the two pericentres, at ~ 0.60 and ~ 1.70 Gyr after the beginning of the 
simulation. The final system is also viewed "edge-on" (bottom-right), showing an axial 
ratio of roughly 1:2. The "young" stars have a somewhat steeper density profile and a 
slightly higher ellipticity. 

elliptical galaxies with higher velocity dispersions than the ones discussed in the Letter. Future 
work may try to understand how such collisions are associated with bigger ellipticals preferentially 
in group centers (e.g. due to two-body relaxation and orbit circularization by dynamical friction). 

How sensitive are our results to the choice of parabolic orbits in the simulations? We extracted 
our typical orbits from cosmological simulations.^^ We have performed a few tests of different 
orbits and orientations, but have not obtained conclusive results about the trends. What we can 
say at this point is that the sample of simulations that we have used spans a range of properties 
that matches many of the real galaxy properties, and especially reproduces many cases of low cjp. 
The dependence on orbit is worth further study. 

2 Robustness to dissipation 

We find that the role of dissipation on the orbit anisotropy beyond i?efr is negligible. The main 
result of dissipation is a more centrally concentrated stellar distribution, which results in ~ 10% 
reduction in velocity dispersion - a minor effect. 
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Figure 4. The tidal origin of the stellar halo and radial orbits is demonstrated via Snap- 
shots from the fiducial Sbc simulation (same as Fig. 3). Marked red is a sample of ("old") 
stars identified at the final time (t = 3) as loosely bound. These stars are shown in their 
positions at several stages along the merger process. The final i?eff is about two small tick 
marks. At the initial time, these stars are spread throughout the progenitor disks {t = 0). 
The first pericentre passage is at t ~ 0.6, and during the rebound (t = 1.0) these stars 
are divided into two groups: those that are flying outward with the extended tidal tails, 
and those that are still associated with the progenitor cores as well as the tidal bridge 
between them. At the second-pericentre coalescence {t ~ 1.8), the outer tidal tails are not 
seen while the inner population is expanding from the center into the halo on highly radial 
orbits. Att = 2.2 we see the stars from the tidal tails falling back in on highly radial orbits, 
to eventually join the inner group in a loosely bound stellar halo extending to a few i?eff- 

One demonstration is provided by the comparison of the "young" stars, involving dissipation, 
to the dissipationless "old" stars. The /? profiles are similar (Fig. 1 of the Letter) and the ap 
profiles differ by only ~ 9% at 3i?efr (Fig. 2 of the Letter). 

A further demonstration is provided in Fig. 5, which compares the results for three merger 
simulations with the initial gas/baryon fraction ranging from to 70%. In particular, we added 
to our original sample a case identical to our fiducial Sbc merger (52% gas) but with all the gas 
particles treated as non-dissipative particles. This "dry" case ends up with a slightly higher (3 
(!), which is balanced by a slightly lower a + 7, resulting in a slightly higher Up. 

We conclude that the radial orbits are generic to mergers independent of whether they are 
"wet" or "dry". The effect of dissipation on Cp, via its effect on the density profile, is limited to 
~ 10%, and is therefore of only secondary concern. This indicates that the key results from our 
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Figure 5. Robustness to dissipation. Anisotropy and velocity dispersion of the "old" stars 
for mergers with different gas fractions. Three projections at two time steps are stacked 
at Res- The red and brown correspond to the same Sbc merger with initially 52% and 
zero gas/baryons respectively. The blue is an Sc merger with 70% initial gas. The effect of 
dissipation is small, at the level of ~ 10% in cXp. 



merger sample are representative of a variety of mergers, "wet" and "dry" , expected in a LCDM 
universe. 

One could argue in passing that the higher a + 7 resulting from the "wet" mergers, being 
closer to the observed a + 7 in the 4 ellipticals discussed (based on Fig. 2 of the Letter), indicates 
that some degree of dissipation is actually required in order to produce the correct density and 
velocity profiles of elliptical galaxies. This is, however, only a marginal indication which is worth 
further study. 

The effects of varying the recipes for star formation and feedback on the relevant features 
addressed in this Letter are yet to be explored. 



3 Robustness to major versus minor mergers 



Our key results are not restricted to 1:1 mergers either. This is demonstrated in Fig. 6, which 
compares our fiducial 1:1 case with several 1:3 and 1:10 mergers. After a single minor merger, the 
remnant does not necessarily resemble an elliptical galaxy (e.g. the disk of the big progenitor is 
still intact), but the halo stars originating in the small progenitor are on very radial orbits, yielding 
final dispersion profiles which are hardly distinguishable from the 1:1 case. The radial orbits and 
the low velocity dispersion are thus generic to the stellar haloes of merger remnants, major or 
minor. This confirms what we heard privately from J. Navarro regarding their preliminary findings 
based on a cosmological sequence of mergers. It indicates that the key results from our merger 
sample are representative of a variety of mergers, major and minor, expected in a LCDM universe. 

Our methodology for checking the effects of minor mergers is limited by the fact that we have 
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Figure 6. Robustness to major versus minor mergers. Anisotropy and velocity dispersion 
for major and minor mergers with different mass ratios. For each merger, three projections 
at two time steps are stacked at Res- Red: our fiducial 1:1 Sbc merger. Blue: four 1:3 G 
mergers stacked together. Green: four 1:10 minor mergers stacked together. In the unequal 
mass mergers, only the "old" stars of the smaller progenitors are used, to better mimic the 
halo population after a sequence of such mergers and avoid the surviving disk component 
of the bigger progenitor. The orbits tend to be more radial for the unequal mass mergers, 
while the final dispersion profiles are hardly distinguishable from the 1:1 case. 

performed so far single mergers only. Our logic is that after multiple minor mergers, the final 
profile is dominated by the many small progenitors which are all contributing to increased radial 
anisotropy. The detailed study of minor mergers is beyond the scope of this Letter. 



4 Robustness to the presence of a bulge 

The effect of a relatively small bulge is seen in Fig. 7. It reduces f3, but increases the central 
stellar concentration, together leaving the velocity dispersion hardly changed. 

We cannot tell from our current simulations what happens when two ellipticals merge. We 
may speculate that if (3 turns out to decrease in such mergers (as indicated by the effect of a bulge), 
while the density profile does not steepen (unlike the case of a bulge), the result may be a slight 
increase in dp. Given that galaxies in groups are more likely to be ellipticals than in the field, 
the dominant group elliptical is more likely to be formed by the mergers of elliptical galaxies (or 
bulge-rich spirals). If /3 is indeed lower in the corresponding remannts, and if the density profiles 
are the same as in spiral mergers, one would predict that dominant group members, which are 
also the most massive ellipticals, should have higher dispersions, as is indirectly observed. This 
will be interesting to address in future work, but it is not required for the main point of our 
current Letter concerning the puzzling low-dispersions in intermediate (L*) ellipticals. 
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Figure 7. Robustness to the presence of a bulge. Anisotropy and velocity dispersion of 
the "old" stars for mergers with and without a 22% bulge component. Three projections 
at two time steps are stacked at Rcs- The purple is a G merger with 20% initial gas to 
stars+gas and 22% bulge to stellar disk ratio. The green is the same but excluding the 
bulge. 

5 Long-term stability 

We have continued one of the G mergers for a longer time, for a total of 6Gyr, where the merger 
occurred at ~ 2.5 Gyr. Shown in Fig. 8 is the time evolution of the profiles of /3 and cjp. During 
the first 0.75 Gyr after the merger (3 dark curves), the system is not yet fully relaxed, though 
the high /3 and low dp are already established. During the following 3 Gyr (blue and light-blue 
curves) the profiles are very stable, at least until 3.5 Gyr after the merger. The slight evolution 
is limited to central regions where remnant gas is slowly forming more stars. 

One can evaluate to what extent the simulated merger remnants obey the Jeans equation, 
given their deviations from spherical symmetry, variations in time, and possible other deviations 
from full equilibrium. Given in each remnant our measured values for a{r), P{r), 7(r) and cTr(r'), 
we have used the spherical Jeans equation to estimate the total mass Mjeans('') within radius r, 
in comparison with the true M(r). We find near 3i2eff for the 10 remnants at about 1 Gyr after 
the merger 7Vfjeans('")/-^('") — 0.9 i 0.15. This is one possible measure of the degree of validity of 
the spherical Jeans equation in these systems. 

6 Parametrizing the velocity-dispersion profiles 

One way to quantify the declining velocity dispersion profile is by fitting a power-law with slope 

7 to o"p(rp) in a given radius range. The slope has been determined for each of the 60 individual 
profiles by linear regression on equally-log-spaced interpolation points in the range 0.5-6i?cff- The 
mean and standard deviation are 7 = 0.53 it 0.16 for "all" the stars, and 7 = 0.61 it 0.22 for 
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Figure 8. Long-term stability. Tlie evolution of the /3 and Up profiles in time after the 
merger (t = 2.5 Gyr) till t = 6 Gyr. The curves are plotted at intervals of 0.25 Gyr; they 
start black (dark blue) and become light blue progressively in time. The profiles are stable 
at least until 3.5 Gyr after the merger. 

the "young" stars. The PN data for the 4 observed galaxies, in the same radius range, give 
7p = 0.59 lb 0.08 lb 0.10, for the mean, error on the mean (using error propagation), and scatter. 
There is an agreement between the simulations and the data, both for "all" the stars and for the 
"young" stars. 

The dp profiles of the four observed galaxies were put together in Fig. 2 of the paper for 
illustrative purposes. The scaling procedure, explained in Methods, is robust to the exact method 
of scaling and the range of data used for it. We tried different minimum radii for the stellar data 
used compared to the fiducial 0.2i?efr- The resulting multiplication factors to get dp = 1 at Rcs for 
the 4 galaxies (NGC 4494, 3379, 821, 4697 respectively) are as follows. For r > 0.2Res: (0.00779, 
0.00761, 0.00637, 0.00714), for r > 0.6i?eff: (0.00758, 0.00740, 0.00656, 0.00708), for r > 1.0i?cfr: 
(0.00788, 0.00806, 0.00623, 0.00685). Using the PNs only we obtain: (0.00804, 0.00779, 0.00656, 
0.00681). The relative changes are of a few percent only. 

7 The full LOSVD and the hi moment 

We find that our sample of merger remnants spans a range of LOSVD's similar to the observed 
sample, as traced both by stars and by PNs. This is illustrated first in Fig. 9 which shows the 
LOSVD of three of our merger remnants in comparison with three observed galaxies, including 
negative and positive kurtosis (or /i4). 

We computed the central values of the fourth moment, h^, for the stars in the simulated 
remnants using narrow radial bins, trying to mimic the typical observer's procedure. We made 
sure that the result within R^g is insensitive to the actual bin size (differences less than 10%) once 
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Figure 9. Line-of-sight velocity distribution P(|fp|/o"p) (LOSVD) of simulation versus 
data. The black curve marks a Gaussian distribution. The LOSVD is computed in two 
radial bins, as marked, within which the variation of cxp is limited. The data is for N3379, 
N4697 and N5128 from left to right respectively. The deviations from Gaussian are small, 
with a tendency for negative and positive kurtosis in N4697 and N5128 respectively. The 
curves are LOSVD from three simulations. Left: the fiducial Sbc merger, stacked views 
from three orthogonal directions at the final time, with a near Gaussian LOSVD. Middle: 
a face-on view of a retrograde merger remnant (G3r) showing a negative kurtosis. Right: 
a face on view of another Sbc merger remnant showing a positive kurtosis. We learn that 
a high anisotropy, (3 ~ 0.5 beyond Res, does not necessarily imply strong deviations from 
a Gaussian distribution. 



the bins are significantly narrower than i?eff- Fig. 10 shows the profiles of /14 for all the simulated 
remnants stacked together. The most relevant curve refers to the average over many random lines 
of sights. The central value is /14 ~ 0.03 =b 0.05. The average within R^s is /i4 ~ 0.02 ± 0.04. In 
the outer range relevant to PNs, 1 — 4i?eff , it is /i4 ~ 0.01 ± 0.05. However, more extreme values 
can be obtained if the line of sight happens to be along one of the principal axes. Along the major 
axis the central /i4 could be 0.1 or larger, and along the minor axis it could be ~ —0.7 at some 
radii and even more negative. 

Comparing the gas-poor G mergers (~ 20% gas) and the gas-rich Sbc mergers (~ 50% gas) 
we find similar results for the average /i4 along random lines of sight, but the G cases tend to 
show less extreme values when viewed along the major or minor axes, reflecting less important 
central disk components in these expected. 

These results are consistent with those of other simulations of "dry" mergers, ^^'^^ /14 = 0.02 ± 
0.04. 

The simulated /14 is consistent with the observed /14 for the integrated stellar light in ellipticals. 
This includes the typical values^^ ~ 0.00 =b 0.04 and those for the four galaxies of the current 
Letter.39-41 

Fig. 10 also shows the /i4 profile for NGC 3379 as measured by four different observations. 
Within 0.3-Refr they span the range ~ 0.0 ± 0.05. The SAURON results"^^ are slightly on the 
high side, ranging from ~ 0.04 to ~ 0.08. At r > 0.3i?efr the /14 measurements are noisier; they 
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Figure 10. The 4th moment of the LOSVD. Left: Profiles of for all simulated remnants 
stacked together, "all" and "young" stars. The hne of sight is either random (blue), or 
along one of the principal axes of the inertia tensor (as labeled). The error bars mark la 
scatter. Right: Profiles of for NGC 3379 from 4 different sources. The la range for 
the "random" view of the simulated remnants is shaded yellow. There is good agreement 
between the simulations and the observed galaxies. 
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extend to values above +0.1 and below -0.1. In SAURON, /i4 is 0.080 — 0.085 when averaged 
over circular bins, with plenty of azimuthal scatter which may be due poor spectral resolution. 
The SAURON /14 profile of NGC 821 is similar, with /14 ~ 0.05 within i?efr- The measured /14 
profiles in these two galaxies are thus consistent with the distribution of /i4 profiles in the merger 
remnants. 

Despite the general agreement with the observed values of /i4, we stress that our sample of 
merger simulations is not supposed to be a fair sample of all the mergers that lead to today's 
ellipticals. Clearly, not all ellipticals, maybe not even most of them, had recent wet major mergers 
of disk systems. However, we demonstrated that the key properties relevant to the present Letter, 
such as the steep density profiles, high /3 and low ap at large radii, are robust to wet and dry 
mergers, major and minor. The simulations do not need to be as representative in terms of the 
inner structure, where the wet mergers tend to lead to a density peak and a small disk, as long as 
these inner features do not affect the key features at large radii. In turn, we do not use here the 
simulated properties in the inner regions as a diagnostic for the physics involved in the formation 
of elliptical galaxies, being dry or wet, etc. 

One lesson from our simulations is that a high anisotropy, f3 ~ 0.5 beyond Rqs, does not nec- 
essarily imply strong deviations from a Gaussian distribution. This suggests that the assumed^^ 
tight correlation between /i4 and /? is model dependent and uncertain in several ways: (a) it is 
based on the assumptions of spherical symmetry and a specific density profile, (b) it is based on 
uncertain and limited orbit-library modeling, (c) it assumes that inner stars and outer PNs are 
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Figure 11. The global structure properties of the simulated merger remnants (symbols, 
colored by the merger type) in comparison with the Fundamental Plane distribution of 
elliptical galaxies in SDSS (la and 2a contours). The dispersion velocity a is measured at 
the central regions. The luminosity is derived from the stellar mass assuming an effective 
M/L = 3. 

the same population, and (d) /14 is also related to rotation. We therefore do not draw conclusions 
from the predictions of P based on /14. 

8 Matching other properties of ellipticals 

Fundamental Plane. The merger remnants qualitatively agree with the Fundamental Plane of 
elliptical galaxies. This is demonstrated in Fig. 11. The deviations are typically at the la level. 
This is work in preparation. 

Surface-brightness profile. The simulated stellar surface-brightness profiles resemble the de 
Vaucouleurs profile remarkably well from below i?efj to several i?efr (Top panel of Fig. 2 of the 
Letter). Some of our simulations produce a central stellar excess (< Ikpc) due to gas settling in 
an inner disk at the remnant center (which might be prevented once AGNs are included). This 
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explains the negative values of /3 at small radii in many of the simulations. Such an excess is not 
seen in typical ellipticals, but is seen in observed mergers.^'^^ This inner excess does not seem to 
affect the kinematics of the relevant, loosely bound stars beyond i?cfr- The line-of-sight velocity 
dispersions are not directly affected by the properties of the material at smaller radii. Indeed, 
we find no correlation between the inner excess and the anisotropy of the halo orbits. The small 
excess is more pronounced in the wet Sbc mergers than in the dryer G mergers, and more so for 
the "young" stars than the "old" ones. The fact that the low Up and steep 7 arc also seen for 
"all" the stars and for the G mergers indicates that the central excess has no important effect on 
our results. 

EUipticity. The remnants of the 1:1 mergers provide projected axial ratios near R^fi that typically 
range from 1:1 to 1:2, not unlike observed ellipticals. In 3D, the median configuration is 1:0.95:0.5, 
namely oblate, but there are triaxial and prolate cases as well, in the ballpark of most observational 
findings.^^"^^ Retrograde mergers are naturally somewhat rounder than the more aligned prograde 
mergers. 

Rotation. The simulated remnants may have a rotating gas-dominated disk at their centers, 
but they are typically far from rotational support at i?eff and beyond. For example, our fiducial 

Sbc case has near R^fi an axial ratio c/a ~ 0.5 (or e ~ 0.5). The cdgc-on view gives a maximum 
of V/ao ~ 0.25 at 1.5ReS, well within the range of the estimates for the 4 observed ellipticals: 
V/<^o = 0.15 — 0.29 (A. Romanowsky, private communication). An isotropic case would have 
given V/ao = [e/{l - e)]^/^ ~ 1- Then (V/ao)* = (F/(7o)/(F/ao,isotropic) ~ 0.25, namely, the 
system ellipticity is dominated by the triaxiality of the velocity dispersion tensor, not by rotation 
(consistent with non-dissipative merger simulations^*^). 

In another merger simulation^^ the remnant looks like a spiral almost immediately after the 
merger. This is perhaps due to their somewhat pathological orbit, where the spins and angular 
momenta are aligned. It is true that our Sbc merger remnants have certain disk-like qualities 
because they started out with a high gas fraction and an extended gas disk. However, the G 
remnants have less disks in them, implying that the high P and low (Xp are not a result of the 
diskiness of the remnant. 

Yet another set of merger simulations^^ indicates a tendency for more tangential orbits at 
large radii. This is likely to be due to their wide angle and slow orbits, in which the tidal effects 
at pericenter are weaker than in our more radial and faster collisions. 

Dark-matter density profile. The simulated DM profiles in the relevant range between i?eff 
and a few i?efF are typically of log slope slightly flatter than -2. This is in the ballpark of what 
might be expected from CDM haloes below ~ 0.1i?vir subject to the baryonic dissipative effects, 
as mentioned in the text. (This could be compared to the unperturbed NFW profile of haloes in 
DM-only simulations, which fiattens to ~ —1.5 near ReS-) 

9 Scatter induced by tricLxiality 

Fig. 3 illustrates the triaxial shape of our typical merger remnant. Fig. 12 demonstrates the effects 
of triaxiality on the velocity dispersion profile. 
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Figure 12. Effect of triaxiality and rotation. Projected velocity dispersion in tlie Sbc 
merger remnant of Fig. 3 as viewed "face-on" (lower) and "edge-on" (upper). Notation as 
in Fig. 2 of the Letter. The triaxiality in the velocity dispersion tensor leads to a ~ 20% 
difference between the two projection directions. The o"p profiles from both views of "all" 
the stars are consistent with the observed cXp. The edge-on view of the "young" stars is also 
in good agreement with the observed cTp, while the face-on cXp generally lies below the data. 
While rotation is part of the anisotropy of the tangential component, the simulated system 
is not supported by rotation: its edge-on rotation velocity is only ~ 25% of the circular 
velocity, with V/ctq ~ 0.2 at ^ _Rofr, and with y/cr ~ 0.3 below Res rising to V/cr 0.5 
beyond 2Res- 

While the four ellipticals addressed here have low, declining velocity dispersion profiles, there 
are other ellipticals which show higher or flatter dispersion profiles. As mention in the text, such 
profiles may be the result of: (a) a specific projection angle in a triaxial or rotating system, (b) 
a high angular-momentum merger, (c) a retrograde merger, (d) a merger with less dissipation, or 
(e) probing with tracers different from PNs, such as globular clusters (of a flatter density profile 
and lower anisotropy). The high Up cases may be associated with more massive galaxies and with 
centers of groups — a subject of future work. 

10 On the PN age 

One may be concerned about the validity of the assertion that the bright PNs are likely to be 
younger than 3 Gyr based on the theoretical work of Marigo et al..^'' 

We should stress first that the issue of whether the PNs are young or not is secondary to our 
results, corresponding to only a small additional reduction of less than 10% in Cp. This is smaller 
than the uncertainties and variation among the galaxies, and smaller than the variations due to 
the merger-orbit parameters or projection direction. Our key result, of radial anisotropy, is even 
more robust, as seen in Fig. 1 of the Letter and as demonstrated above. 
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Still, the theoretical analysis seems straightforward and we are not aware of convincing theo- 
retical arguments against its validity. The similar luminosity function of PNs in different ellipticals 
is an empirical result which should rather be explained (e.g. by a similar merger history for these 
galaxies, including many minor mergers), but it docs not seem to justify ignoring the theoret- 
ical predictions based on stellar theory. It is therefore fair to consider the PN young age as a 
possibility. 

We may comment in passing that the hypothesis that the PNs are blue stragglers produced 
by stellar collisions, may perhaps work for globular clusters, but it is hard to see how the stellar 
collision rate could be high enough at the outskirts of elliptical galaxies. Perhaps a recent merger 
would help... The idea of blue stragglers may work if they are produced by mass transfer in 
common binary systems. 

11 Rising versus constant P 

Is the gradual rise in /3(r) responsible for a steeper decline in cjp? This seems consistent with the 
trend seen in the comparison of "dry" and "wet" simulations in Fig. 5, for which /? is a constant 
and rising respectively. 

However, a test using the spherical Jeans equation indicates that low /3 values below i?cfr have 
negligible effects on c7p beyond Res- Fig. 13 compares two solutions of the Jeans equation, one 
with P = 0.5 at all radii, and the other with /3 = at small radii, smoothly rising near R^s to 
P = 0.5 at large radii following an assumed functional form.^^ For example, f3 = 0.25 at 1.4ii!eff- 
The shape of the /? profile is similar to the average shape in our simulated remnants. After 
matching the solutions at R^s, the dp profiles are almost identical at all larger radii. Based on 
this example, we suspect that the rising nature of the profile may not to be the key factor driving 
the decline in (jp. 

We recently learned from A. Romanowsky (private comm.) about the /3 profile in their 
preferred orbit-library solution to N3379: it has a negative value of /3 = —0.3 at -Rcflf; and it 
rises toward ~ -1-0.5 but only beyond 5i?eff- While being rising in part of the relevant range, it 
does not seem to reach high enough values there, and therefore does not support a massive DM 
halo. In comparison, in many of our simulations, /? is positive at ReS, it is close to 0.5 already at 
2i?eff, and it keeps rising toward larger radii. 

12 Why was the DM solution missed before? 

Using merger simulations, we provide "a" fully specified solution, consistent with what we know 
about DM haloes in other galaxies and in the ACDM scenario, which seems not to violate all the 
relevant observed aspects within the given uncertainties. Independent of whether this solution 
is indeed "the" solution for elliptical galaxies, it is claimed^ (R03) to have been missed in the 
original application of the orbit-library maximum-likelihood algorithm. 

We cannot provide a definite answer to why this solution was missed; the published material 
(as well as private communication with the authors) provided only partial information for evalu- 
ating the method or for fully characterizing their preferred solutions. A sample output from our 
simulations can provide a useful testbed for their method. 
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Figure 13. Rising versus constant f3. A comparison of solutions to the spherical Jeans 
equation; one with (3 = 0.5 at all radii and the other with a rising (3{r) from zero at small 
radii to 0.5 at large radii, with the transition near i?eff- 

However, it is not at all surprising that a complex maximum-likelihood method based on an 
orbit library may miss some of the solutions. We can only raise several ideas for where such 
method could go wrong, as follows: 

1. Assuming a stationary potential. The outer system may still be evolving (see Fig. 4). We 
find in our Sbc simulations that the Jeans equation is typically valid to ~ 10% accuracy, in 
the sense that the total mass within 1 — 3Res as derived by the Jeans equation is ~ 90% of 
the true value. In a few cases the deviation can be larger. 

2. Assuming spherical symmetry. Ellipticals can appear circular in one projection while actu- 
ally having an axial ratio of 1:2 (see Fig. 3). 

3. Assuming density profiles with a specific, limited shape, such as NFW for the DM. There 
is no reason to believe that the DM haloes in ellipticals, under the influence of dissipation 
within them, still obey their dissipationless NFW shape near 

i2efT-^^ In fact, our simula- 
tions demonstrate that a major merger of NFW haloes including disks may lead to slight 
deviations from NFW in the relevant range (see Fig. 1 of the Letter, showing a power-law 
of slope ~ —2 in the relevant range). 

4. Their orbit-library method heavily depends on constraints from the values of higher velocity 
moments, such as /14, as determined from the central stellar diffuse light. These constraints 
may be misleading if the PNs detected at large radii actually represent a population that 
is kinematically different from the central stars. 

5. The problem may lie in some incompleteness of the orbit library used, or in the way by which 
it is composed for best fit. There were claims in the literature about failure to converge 
when the number of orbits is increased. ^^'^^'^^ One may suspect that the R03 method is 
subject to such problems. 




6. 



A multi-parameter (or "non-parametric") maximum-likelihood method may easily miss 
some of the solutions. 
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13 Robustness and uniqueness of the results 

Unfortunately, full exact predictions of the ACDM scenario are beyond what can be achieved 
by current methods. We therefore use state-of-the-art simulations to study in detail a sample 
of merger cases. This sample, based on plausible assumptions, provides DM-dominated systems 
which are consistent with the observed ellipticals in all relevant testable respects, and it repro- 
duces the puzzling low velocity dispersions and their variations among different galaxies. This by 
itself is an important step forward, given that previous work failed to find such a solution and 
thus challenged the whole cosmological scenario. Demonstrating that the low velocities are not 
necessarily fatal for ACDM is a step forward by itself. 

However, the variety of cases simulated provide reasonably compelling evidence for the robust- 
ness of the key phenomenon — radial orbits — indicating that this is indeed a generic outcome 
of a large class of mergers in CDM-like scenarios. In particular, we have demonstrated robustness 
to dissipation, mass ratio, the presence of a bulge and several other merger characteristics, and 
have gained insight into the generic origin of the effect via tidal processes. 

An alternative solution has been proposed, appealing to a low halo concentration parameter 
C. Using our Jeans-equation algorithm, we find that even a C as small as 30% of the mean 
ACDM value fails to fit the ap profile of NGC 3379 unless (3 ^ 0.5 near ReS and beyond. The 
fit proposed^^ requires an unrealistic C, five times smaller than the ACDM value. This indicates 
that our solution, based on radial stellar orbits, is the only viable solution proposed so far in the 
literature within the ACDM cosmology. 

We have demonstrated that the solutions emerging from our sample merger simulations are 
representative of a large class of mergers likely to occur in a ACDM cosmology. We conclude 
that the radial anisotropics necessary for compatibility with the low ap observed are not really 
"pathological" but generally expected. Then, triaxiality of the galaxies could play a major role 
in observing low dispersions; recent star formation could in rare cases also contribute." 
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